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Abstract 

The Strecker reaction of acetaldehyde, NH3, and HCN to afford alanine was studied by DFT calculations for the first time, which 
involves two reaction stages. In the first reaction stage, the aminonitrile was formed. The rate-determining step is the deprotonation 
of the NH3" 1 " group in MeCH(OH)-NH3 + to form 1-aminoethanol, which occurs with an activation energy barrier (AE^) of 9.6 kcal/ 
mol. The stereochemistry (R or 5) of the aminonitrile product is determined at the NH3 addition step to the carbonyl carbon of the 
aldehyde. While the addition of CN~ to the carbon atom of the protonated imine 7 appears to scramble the stereochemistry, the 
water cluster above the imine plane reinforces the CN~ to attack the imine group below the plane. The enforcement hinders the 
scrambling. In the second stage, the aminonitrile transforms to alanine, where an amide Me-CH(NH2)-C(=0)-NH2 is the key inter- 
mediate. The rate-determining step is the hydrolysis of the cyano group of N(amino)-protonated aminonitrile which occurs with an 
AE^ value of 34.7 kcal/mol. In the Strecker reaction, the proton transfer along the hydrogen bonds plays a crucial role. 



Introduction 

In 1850, Adolph Strecker reported a reaction that affords 
alanine from acetaldehyde, ammonia and hydrogen cyanide [1]. 
The original form of Strecker amino acid synthesis is shown in 
Scheme 1(a). In this reaction, the aldehyde reacts with hydrogen 
cyanide to form an aminonitrile, which undergoes hydrolysis to 
afford alanine in the acidic solution. The traditional Strecker 
reaction gave racemic a-aminonitriles (mixtures of equal 
amounts of R and S forms), where an imine RCH=NH was 



considered to be the key intermediate [2]. Three typical reac- 
tions are presented in Scheme 1(b) [3]. 

In 1963, Harada reported the first asymmetric Strecker reaction, 
in which an (S)-a-phenylethylamine was employed as the chiral 
auxiliary [4]. In this reaction, he obtained a chiral alanine with 
95% optically activity; see Scheme 2. In 1996, Lipton et al. 
succeeded in a series of asymmetric Strecker reactions by 
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Scheme 1 : The general form of the Strecker reaction. The reaction (b) is taken from [2]. 
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employing a chiral catalyst, a cyclic dipeptide [5]. In these reac- However, when benzaldehyde and NH 3 instead of the N-substi- 

tions, N-substituted imines react with HCN to yield (S)-a- tuted imine were employed as the substrates, the reaction 

aminonitriles with remarkably high enantiomeric excess (ee). afforded an initial product Ph-CH(NH2)-CN of configurational 

One example is shown in Scheme 3. instability[5]. In the following, Sigman and Jacobsen used a 
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Scheme 3: The first asymmetric synthesis of a-aminonitirles via a chiral catalyst [5]. 
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parallel combinatorial library synthesis for the discovery and 
optimization of a chiral catalyst for the reaction of imines and 
HCN [6]. From then on, various catalytic asymmetric Strecker 
reactions have been reported to gain high enantioselectivity of 
the hydrocyanation reaction of imines [7-12], However, the 
origin of the enantioselectivity in the asymmetric Strecker reac- 
tions has not been clarified. 

To our knowledge, the elementary processes of the whole 
Strecker reaction have not been elucidated. As shown in 
Scheme 1, the Strecker reaction includes two reaction stages. 
The first reaction stage is the condensation of aldehydes with 
ammonia and hydrogen cyanide leading to a-aminonitriles . The 
second reaction stage is the hydrolysis of the nitrile group. In 
these reactions, K + (or Na + ) and Cl~ ions are not involved, as 
shown in Scheme 1(a). Therefore, it is suitable to the theoreti- 
cal investigation of the reaction mechanism, because the effect 
of counter ions does not need to be considered. 

Actually, several theoretical studies were reported of the last 
step of the first reaction stage of the Strecker reaction [10-16], 
i.e. the hydrocyanation of imines (or protonated imines + CN~) 
to aminonitriles. In those works, how the nucleophile CN~ is 
generated has not been examined. Because HCN is a very weak 
acid with a dissociation constant of K a = 1.3 x 10~ 9 mol/L (in 
water, 18 °C), direct dissociation reaction HCN — > H + + CN~ is 
difficult to occur. 

In the second reaction stage, i.e., the acid-catalyzed hydrolysis 
of the cyano group, protonation of the group appears to cause 
the addition of OH 2 to the cyano carbon: 

R-CH(NH 2 )-CN + H + R-CH(NH 2 )-CNH + 

R-CH(NH 2 )-CNH + + OH 2 -> R-CH(NH 2 )-C(OH)=NH + H + 

However, the proton affinity (PA) of the nitrile is much smaller 
than that of the amino group, for example, the PAs of the cyano 
and amino groups of 2-amino-propanonitrile (Me-CH(NH 2 )- 
CN) are 190.7 and 199.6 kcal/mol, respectively. Thus, in the 
acidic solution (2H 2 0 + 2HC1), Me-CH(NH 3 + )-CN should be 
afforded; see Scheme 1(b). The reaction mechanism of this hy- 
drolysis is also unclear. 

To address the above issues, we performed DFT calculations of 
the Strecker reaction shown in Scheme 1(a). Here, ten specific 
water molecules were considered. 

Methods of calculation 

Geometry optimizations were performed by density functional 
theory (DFT) with the B3LYP [17,18] functional. The basis set 



6-311+G(d,p) was employed for all the atoms in the calcula- 
tions. The solution (water) effect was considered by the Polariz- 
able Continuum Model (PCM) [19-21]. Vibrational analyses 
were carried out to make sure whether a stationary point is an 
equilibrium structure or a transition state (TS). From TSs, reac- 
tion paths were traced by the intrinsic reaction coordinate (IRC) 
method [22,23] to obtain the energy-minimum geometries. All 
the calculations were carried out using the GAUSSIAN 09 [24] 
program package. Throughout this paper, the discussion was 
presented based on the potential energy changes with zero-point 
vibrational energy (ZPE) correction unless otherwise noted. 

Results and Discussion 

Formation reaction of aminonitrile (the first 

stage) 

The reaction model is shown in Scheme 4. In the model, lone- 
pair electrons of the oxygen and nitrogen atoms participate in 
hydrogen bonds. In calculating each TS, ten water molecules 
were placed so that the large hydrogen-bond stabilization is 
gained. After geometries of TSs were determined, those of 
energy minima were obtained by IRC and the subsequent opti- 
mizations. By the use of the similar water cluster models, TS 
geometries and activation energies (E^s) in the base promoted 
ester hydrolyses were calculated [25]. In Ph-COOE t + 
OH"(H 2 0)„->Ph-COO" + HO-E t + (H 2 0)„, E a = +14.7 kcal/ 
mol (n = 5), +16.3 {n = 8), +16.3 {n = 12) and +15.6 (n = 16) 
were obtained, where the experimental E & is + 14.6 kcal/mol. 
Also, in the Bamberger rearrangement Ph-NH(OH) + 
(H 3 0 + ) 2 (H 2 0) 13 -y para-HO-C 6 H 4 -NH 3 + + H 3 0 + (H 2 0) 14 [26], 
E a = +26.3 kcal/mol was calculated, where the experimental E a 
is +24.8 kcal/mol. 




Scheme 4: A reaction model composed of Me-CH=0, HCN, NH3 and 
(H20)io for geometry optimizations to trace elementary processes. 
Broken lines stand for hydrogen bonds. 
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Scheme 5: Possible pathways for the formation of aminonitrile from acetaldehyde. 



Our proposed reaction pathways are shown in Scheme 5. 
Geometries of TSs are shown in Figure 1 and those of precursor 
1, intermediates and product 8 are provided in Supporting Infor- 
mation File 1 Figure SI. 

As shown in Figure SI, Supporting Information File 1, 
MeCH=0, NH 3 and HCN are separated by the H 2 0 cluster in 
the precursor complex 1. The reaction begins with the addition 
of NH3 to the carbonyl carbon of acetaldehyde to form a 
Mulliken charge-transfer complex 2. This complex was firstly 
proposed here. In 2, the C-0 bond is elongated to 1.345 A, 
which has an alkoxide character and the complex is not stable in 
the gas phase. However, it is more stable than the precursor 
complex by 3.5 kcal/mol when ten specific water molecules are 
considered; see Figure 2. This result indicates that the consider- 
ation of water molecules in the reaction is necessary to describe 
the step 1 — > 2. Then, the alkoxide oxygen atom captures a 
proton from a surrounding water molecule to form 
Me(H)C(OH)-NH 3 + and a remaining OH" anion in the 
surrounding. They form an ion pair 3 [Me(H)C(OH)-NH3 + and 
OH - ]. Next, the hydroxide ion catches a proton from HCN 
through the transition state TS3/4 to afford a more stable ion- 
pair intermediate 4 [Me(H)C(OH)-NH 3 + and CN"]. This step is 
exothermic by 5.5 kcal/mol with a small activation energy 
barrier (AE^) of 1.4 kcal/mol. Starting from 4, a proton of the 
NH 3 + group migrates to one water molecule to form 5 
[Me(H)C(OH)-NH 2 (H 3 0 + ) and CN"] via a transition state 
TS4/5. After that, the proton migrates from H 3 0 + to the hydroxy 
group of Me(H)C(OH)-NH 2 via TS 5/6 to yield Me(H)C(OH 2 + )- 
NH 2 6. From 6, H 2 0 is easily eliminated through TSj/7 to 
afford the protonated imine 7 with a AE^ value of 0.3 kcal/mol. 
At last, CN~ nucleophilically attacks the carbon atom of 



MeCH=NH 2 through TS7/8 to afford a 2-aminopropanonitrile 
8. 

Starting from 4, a concerted S^-type pathway was also exam- 
ined, which directly leads to the nitrile compound 8; see 
Scheme 6. In this pathway, the proton is transferred from the 
NH 3 + group to the hydroxy group via a two-water-molecule 
bridge. At the same time, the H 2 0 elimination and the ap- 
proach of CN~ to 4 concomitantly take place with a Walden 
inversion. The transition state TS4/8 was successfully located; 
see Figure 1. However, this pathway needs a large energy 
barrier of 58.8 = [+51.3 - (-7.5)] kcal/mol, indicating that it is 
difficult to occur. 

Starting from the ion-pair intermediate 3, we also investigated 
the possibility that the hydroxide ion captures a proton from the 
NH 3 + group to form 9, Me-C(H)(OH)-NH 3 + + OH" -> 
Me-C(H)(OH)-NH 2 + OH 2 . This reaction step occurs through a 
transition state TS3/9 with a A£^ value of 1 .9 kcal/mol, which is 
comparable to that (1.4 kcal/mol) of the proton transfer step 




Scheme 6: A short-cut path by the nucleophilic displacement and the 
concomitant proton transfer. "The first bypass" in Scheme 5. 
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Figure 1: Geometries of transition states along the reaction from acetaldehyde (1) to the aminonitrile 8. Distances are in A. TS1/2 means, for instance, 
a transition state for the step 1 — > 2. 



from HCN to OH~. However, 9 is less stable than 4 by 1.8 kcal/ 
mol, indicating that the OH~ prefers to capture a proton from 
HCN rather than from the NH3 + group. 

As shown in Figure 2, the rate-determining step of this reaction 
stage is the proton migration from the NHj + group to the water 
cluster (from 4 to 5), where the energy barrier is 17.1= [+9.6 - 
(-7.5)] kcal/mol. Other proton transfer steps facilely occur. 
These results are in consistent with the room-temperature 
experimental condition in Scheme 1 . For the rate-determining 
step, we also checked an extended model "TS4/5-ext", where 
ten water molecules are added (the molecular formula, 
C3H48N2O21). The geometry of TS4/5-ext is shown in 



Supporting Information File 1 Figure S2. The geometrical para- 
meters of the proton-transfer region of TS4/5-ext are similar to 
that of TS4/5. Also, the energy difference between TS4/5-ext 
and 4-ext (17.3 kcal/mol) is very close to that (17.1 kcal/mol) 
between TS4/5 and 4. 

As shown in Figure 1, the product, aminonitrile 8, is in an S 
form. However, racemic a-aminonitriles are obtained experi- 
mentally. This stereochemical scrambling is explicable on the 
basis of the computational results, as follows. In TS1/2, the NH3 
molecule may add to MeCH=0 from both upper and lower 
directions equivalently, which leads to the racemic products. 
However, in TS7/8, the nucleophile CN~ is obligated to attack 
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Figure 2: Energy changes along elementary processes from acetaldehyde to aminonitrile. Bold numbers are defined in Scheme 5. 



MeC(H)=NH2 + at the plane opposite to the OH2 dissociating 
side (see Scheme 7). The addition model of TS1/2 was exam- 
ined by the use of the amine in Scheme 2. The activation energy 
of the less crowded TS1/2— A ' s 1-8 kcal/mol smaller than the 
more crowded TS1/2-B; see Figure 3. This calculation result is 
consistent with Harada's work that a chiral product [4] was 
obtained in the Strecker reaction. 

Hydrolysis of amino nitrile to amino acid (the 
second stage) 

In the acidic hydrolysis of amino nitrile, we take 2-amino- 
propanonitrile +H30 + (H20)io (8) as a precursor complex; see 
Figure S3 in Supporting Information File 1 for the geometry of 
8. Our proposed pathways are shown in Scheme 8. 




Figure 3: Two transition states (A and B) of the nucleophilic addition of (S)-a-phenylethylamine to acetaldehyde. (H20)io is also included, and the 
molecular formula of the reaction system is C11H36N2O11. 
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Scheme 7: A contrast of the nucleophilic addition. 
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There are two competitive pathways from the precursor com- 
plex: One (I) is the protonation of the amino group to form a 
N(on amino)-protonated aminonitrile 10. This step occurs 
through TSg/io with the AE* and AE values of 4.1 and -6.8 
kcal/mol, respectively; see Figure 4. 

The other pathway (II) is the hydrolysis of the C=N group to 
form a compound 18 through TSg/is, where the OH group is 
added to the carbon atom and the hydrogen atom attachs to the 
nitrogen atom. This step needs a considerably large energy 
barrier of 32.0 kcal/mol with an endothermicity of 0.2 kcal/mol. 
Obviously, the protonation of the amino group (I) is much more 
favorable than the hydrolysis of the cyano group (II). As a 
result, compound 10 is the starting point for the following reac- 
tions. From 10, a water trimer reacts with the cyano group 
through TSio/n to afford the N(on amino)-protonated 2-amino- 
1-hydroxypropanimine 11. At TSio/n, the hydroxy group adds 
to the carbon atom of the cyano carbon nucleophilically. Simul- 
taneously, the proton migrates to the nitrogen atom through a 
two-water-bridge; see Figure 5 for the geometry of TSio/n- The 
C=N group in 10 convers to a C(OH)=NH group in 11. This 
reaction step is endothermic by 4.6 = [-2.2 - (-6.8)] kcal/mol 



with a large AE* value of 34.7 = [+27.9 - (-6.8)] kcal/mol. 
Although this AE* value is apparently larger than that of TSs/is, 
TSio/n lies lower than TSs/is by 4.1 kcal/mol when taking the 
energy of the precursor complex as a reference; see Figure 4. 
We examined the role of the NH3 + group in the reaction by 
investigating the hydrolysis of the cyano group of a methyl- 
substituted model 10(Me). In this model, we replaced the NH 3 + 
group in 10 with a methyl group. The AE* value of the hydroly- 
sis step increases to 34.6 kcal/mol. It indicates that the NH3 + 
group enhances the electrophilicity of the cyano carbon, which 
is favorable for the OH2 addition. The TS geometry of the 
methyl substituted model, TSio/n(Me), is shown in Supporting 
Information File 1 Figure S6. 

In the following, a proton on the NH3" 1 " group in 11 is trans- 
ferred to the imine nitrogen to afford 12. This proton transfer 
step is facilitated by a water molecule bridge. The AE^ and AE 
values are 5.5 = [+3.3 - (-2.2)] and -3.8 = [-6.0 -(2.2)] kcal/ 
mol, respectively. Starting from 12, there are two possible path- 
ways, paths A and B, to form the final alanine product; see 
Scheme 8. In path A, the deprotonation of the amino group 
occurs to produce an amide intermediate 13 with an exother- 
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Figure 4: Energy changes along elementary processes from 2-amino nitrile 8 to 2-amino acid 16. Brown-color lines stand for the most favorable 
route. 
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Figure 5: Geometries of transition states along the most favorable route from 2-aminonitrile 8 to 2-amino acid 16. 
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energy of 13 are 30.5 and 25.8 kcal/mol, respectively. After 
that, the NH3 moiety is ready to dissociate from 15 to produce 
the product (i?)-alanine 16 with a AE^ value of 0.7 kcal/mol of 
TS15/16. In path B, the second H2O molecule is added to the 
C=N double bond in 12 to form an intermediate 17 through 
TS12/17. The AE^ and AE values of this step are 25.5 and 8.6 
kcal/mol, respectively. However, TS12/17 li es higher than 
TS14/15 by 3.1 kcal/mol and 17 in path B is much more unstable 
than the intermediate 13 in path A by 17.0 kcal/mol. These 
energy differences suggest that the deprotonation of the amino 
group occurs more favorably. Then, from 17 the elimination of 
the NH4 + group takes place to afford the alanine 16 with a Afi^ 
value of 2.5 kcal/mol. According to the above discussion, the 
path-B is less favorable than the path A. 

The most favorable pathway for the second stage of the Strecker 
reaction was shown in brown color in Figure 4. The rate-deter- 
mining step is the OH 2 addition to the cyano group (TSio/n) 
with the activation energy of 27.9 kcal/mol. The competitive 
TSg/18 has an energy barrier of 32.0 kcal/mol. Geometries along 
the unfavorable routes pathway II and pathway IB, are shown 
in Figures S4 and S5, Supporting Information File 1, respective- 
ly. The relative stability of these two transition states were 
checked with extended models TSjo/n-ext and TSg/ig-ext, 
which have a molecular formula of C3H49N2C>2i + . Their 



geometries are shown in Figure S7. TSg/ig-ext lies higher than 
TSio/n-ext by 3.4 kcal/mol, which is consistent with the small 
model system. 

Conclusion 

In this work, we theoretically investigated the whole Strecker 
reaction shown in Scheme 1(a), which includes two reaction 
stages. The most favorable pathways are summarized in 
Scheme 9. 

As shown in the upper half of Scheme 9, the first reaction stage, 
acetaldehyde + NH 3 + HCN + (H 2 O) 10 (1) -» 2-amino- 
propanonitrile (H 2 0)n(8), is composed of seven elementary 
processes. The rate-determining step is the deprotonation of the 
NH3 + group in MeCH(OH)-NH3 + to form 1-aminoethanol, 
which occurs with an activation energy barrier of 9.6 kcal/mol. 
The stereochemistry (R or S) of the product aminonitrile is 
determined by equal addition of NH 3 to the carbonyl carbon of 
the aldehyde in both sides. While the addition of CN" to the 
carbon atom of the protonated imine 7 appears to give the 
scrambling of the stereochemistry, the water cluster above the 
imine plane reinforces the CN" to attack the carbon atom below 
the plane; see Scheme 7. While HCN is a very weak acid, CN~ 
may be generated by the proton transfer, HCN + OH~ — > CN~ + 
H2O (3 — > TS3/4 — » 4) in this reaction stage. As shown in the 
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The second stage: hydrolysis of aminonitrile to form amino acid 
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Scheme 9: Summary of the present computational work expressed by minimal models. 
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lower half of Scheme 9, the second reaction stage, aminonitrile 
+ H 3 0 + (H 2 0)io — » alanine + NH 4 + (H 2 0)9, is also composed of 
seven elementary processes. In this reaction stage, the protona- 
tion to the amino nitrogen occurs first, which enhances the 
subsequent hydrolysis of the cyano group to form an imine 
C(OH)=NH moiety. Its rate-determining step is the hydrolysis 
of the cyano group of N(amino)-protonated aminonitrile to 
afford a N(amino)-protonated 1 -hydro xypropanimine 11. The 
AE^ value of this step is 34.7 kcal/mol and the large value 
corresponds to the high temperature conditions in Scheme 1(b). 
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